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ABSTRACT 

Nomograms of mean column heating due to surface sensible and latent heat fluxes have 
been developed from Stage and Businger’s (1981a, b) boundary layer model for cold air out- 
breaks over warm water. Mean sensible heating of the cloud free region is related to the cloud 
free path (CFP, the distance fram the shore to the first cloud formation) and the difference be- 
tween land air and sea surface temperatures, 6 j and Oq, respectively. Mean latent heating is re- 
lated to the CFP and the difference between land air and sea surface specific humidities, qj and 
q 0 , respectively. Results are also applicable to any path within the cloud free region. Corres- 
ponding heat fluxes may be obtained by multiplying the mean heating by the mean wind speed 
in the boundary layer, The sensible heating estimated by the present method is found to be in 
good agreement with that computed from the bulk transfer formula. The sensitivity of the solu- 
tions to the variations in the initial coastal soundings and large scale subsidence is also investigated, 
Tiie results are not sensitive to divergence but are affected by the inital lapse rate of potential 
temperature; the greater the stability, the smaller the heating, other things being equal, Unless 
one knows the lapse rate at the shore, this requires another independent measurement. For this 
purpose we propose to use the downwind slope of the square of the boundary layer height, the 
mean value of which is also directly proportional to the mean sensible heating. The height of the 
boundary layer should be measurable by future spacebome lidar systems. The general behavior 
of the mean sensible heating, the potential temperature and the height of the boundary layer as 
a function of downwind distance within the cloud free region and their relations to several impor- 
tant parameters are studied analytically in the Appendix. By-products include the finding that 
the sensible (latent) heat flux is virtually linear with the contrast in land air and sea surface tem- 
perature (specific humidities), thus providing a new kind of flux parameterization in lieu of the 
classical bulk transfer formulas, The applicability of the results to lake-effect snowstorms is also 
noted. Finally, the method can be used in reverse to check the validity of boundary layer models. 
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ESTIMATING OCEAN-AIR HEAT FLUXES 
DURING COLD AIR OUTBREAKS BY SATELLITE 

I. INTRODUCTION 


Airmen interaction plays a very important role in a wide range of weather, climate, and 
ocean problems. There is increasingly persuasive evidence that anomalies in sea surface temper- 
ature (SST) control both regional and global climatic anomalies (e.g,, Bjerknes, 1966, 1969; 
Namias, 1969, 1978), The anomalous SST obviously affects the atmospheric circulation through 
the ocean-air heat fluxes. For this reason much effort is going into the accurate parameterization 
of the fluxes of heat, moisture, and momentum in GCM models. Tire ocean to air heat fluxes 
are also related to the prevailing synoptic pattern and mcsoscatc convective system (Rung and 
Siegal, 1979; Sheu and Agee, 1977; Warsh, 1973). These fluxes over the warm water of the 
Kuroshio current and the Gulf stream are generally known to have profound effects on cyclogen- 
esis and storm intensification, particularly on the development of the wintertime cyclones of the 
Taiwan and Cape Hatteras lows (Agee and Howley, 1977), 


Air-sea heat exchange obviously has significant manifestations within the sea as well. For 
example, Worthington (1977) has postulated that the Gulf stream accelerates after a cold winter 
as a result of the enhanced differential cross-stream cooling and density gradient. On the other 
hand Huh 1 has found that the winter heat loss from the Gulf of Mexico takes place episodically 
during cold air outbreaks and that there is a virtual balance between the cooling of the coastal 
waters and the heat given up to the atmosphere* 

Possible approaches toward the measurement of ocean-air heat fluxes were treated at length 
at the Jet Propulsion Laboratory/Scripps Institute of Oceanography Workshop on air-sea interac- 
tion (NASA, 1980) and the Workshop on Ocean Surface Energetics (Gautier, 1981). One of the 
key conclusions of the former meeting was that the prospects for the measurement of the fluxes 

'Private communication from Professor Oscar Huh, Louisiana State University 
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of heat, moisture, and momentum from space are not promising. This is due to the fact that 
such fluxes depend upon the profiles of wind, temperature, and moisture in the atmospheric 
boundary layer, or equivalently, on the bulk transport equations, and there are no clear methods 
of sensing the required parameters remotely, Only slightly greater optimism was expressed at 
the second workshop mentioned above. 

On the other hand, the fluxes of heat and moisture from sea to air often produce effects 
which can be observed and which do provide some means of estimating these fluxes, The most 
obvious products of such fluxes are clouds, the growth of the convective boundary layer, and 
the warming and moistening of that layer (e,g„ Lenchow, 1973). Such effects are especially dra- 
matic in the case of cold air outbreaks over warm waters along the east coasts of Asia and North 
America (e,g„ Ninomiya, 1976; Ninomiya and Akiyamu, 1976). There the effects arc clearly 
manifested in the form of cloud streets, m shown in Fig, !, We see that the upwind (i.e., near- 
est the shore) edge of the cloud streets is virtually a displaced image of the coastline, the dis- 
placement. being along the wind direction. It appeared likely that this displacement or distance 
between shore and the upwind edge of the cloud, here called the “cloud free path” or CF1\ is 
related to the heat and moisture fluxes from the coastal waters such that the greater the fluxes 
the smaller the CFP. Tins was essentially our starting hypothesis. In this paper we show that 
this is indeed the case and develop a method by which one may estimate the fluxes with the aid 
of the CFP, the sea surface temperature and saturation water vapor mixing ratio, and certain 
other quantities to he described later. 

While the limitations of this study to cold air outbreaks and legions within 200 to 300 km 
of the coasts would appear rather restrictive, these areas and periods include the major share of 
the annual global ocean heat losses (Bunker and Worthington, 1976; Agee and Howley, 1977). 
Kung and Siegel further suggest that the heat transfer from the warm water to the atmosphere 
in these regions plays a distinct role in the Northern Hemisphere winter circulation, 
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2. BACKGROUND 


Henry and Thompson (1976) inadvertently provided a hint us to how one might use satel- 
lite observations to estimate heat fluxes during cold air outbreaks. They used the radiosonde 
temperature and humidity profiles observed at the coast of the Gulf of Mexico near New Orleans 
and 24 hours later at Merida. Yucatan, Mexico to deduce the heat fluxes at various points of the 
trajectory following the cold northerly flow by budget calculation with simple interpolation of 
temperature and humidity soundings. Well developed cloud streets were evident starting about 
100 km south of the Louisiana coast. In their case,, the very cold land air of about -4°C over 
warm Gulf waters of about 20 e C produced a total average heat flux of 2400 Win' 2 in the first 
1 10 km of travel to a point close to the edge of the first clouds. This is a major flux and again 
indicates why the problem is important in the previously noted contexts. We shall show later 
that the fluxes of Henry and Thompson are probably overestimated. 

The question at issue, however, is how one can use satellite observations to estimate the 
heat fluxes without monitoring the temperature and moisture of the air column at both ends of 
the trajectory as was done by Henry and Thompson. In a preliminary version of this paper, At- 
las et al. (1981), hereafter referred to as I, tried several approaches. Their first approach was to 
use the satellite observed cloud top temperature at initial cloud formation. This temperature 
could be extrapolated down to near the surface adiabatically, thus providing an estimated sound- 
ing which could then be differenced from that measured by radiosonde at the shore to give sen- 
sible heating. Similarly, the temperature at cloud formation defines the saturation water vapor 
mixing ratio which is also well mixed down to near the surface, Comparison to values at the 
shore would thus provide latent heating, The key assumption in this approach is that the sound- 
ing along the trajectory remains unchanged above the flux modified convective boundary layer so 
that the satellite measured cloud top temperature could be assigned to a known altitude. 

While the elements of the concept are sound, its implementation is subject to excessive er- 
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rors due largely to the small size of the initial clouds relative to the I km dimension of the 
AVHRR (Advanced Very High Resolution Radiometer) IR Held of view on the TIROS satellites 
and the indeterminacy of the temperature correction for unfilled fields of view, These errors 
are compounded by the large vertical gradient of potential temperature at the inversion so that 
small errors in determining the temperature of the initial cloud tops, assumed to be at the inver- 
sion, translate to large errors in the deduced sounding and thus in net heating. While this method 
does not work with existing observing systems, it should be kept in mind for possible future ap- 
plication with higher resolution IR imagers and lidar sounders capable of independent cloud 
height determination, 

In I we also utilized the results of Stage (1979) who, with the aid of the boundary layer 
model of Stage and Businger (198 la, b), computed the potential temperature, vapor mixing ratio, 
and the height of the convective boundary layer as a function of distance for cold air transports 
over Lake Ontario, lie also calculated the distance from shore to first cloud formation (i.e,, the 
“cloud free path” or CFP), all as a function of varying lake temperatures. It was only necessary 
then to estimate the sensible and latent heating due to surface fluxes and, with the aid of certain 
approximations, develop the desired relations to measurable parameters. Because Stage (1979) 
had only run the Stage-Businger model for a single land air temperature and lapse rate, and for a 
limited range of lake water temperature, it was not possible to generalize the results in 1 to a suf- 
ficiently wide range of conditions or to determine their limitations, In the present paper we re- 
port the generalized results obtained by running the Stage-Businger model over a very broad range 
of conditions and correct some of the erroneous inferences which were made in I on the basis of 
more limited data, 

In particular, we have developed generalized nomograms of mean column heating due to the 
surface sensible and latent heat fluxes. The sensible heat nomogram relates the mean sensible 
heating to the cloud free path (CFP) and the difference between land air and sea surface temper- 
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atures, 0,and 0 O , respectively. The latent heat nomogram relates the mean latent heating to the 
CFP and the difference between land air and sea surface vapor mixing ratio, qj and q 0 , respec- 
tively. The aim, of course, is to estimate the mean sensible and latent heating in the cloud free 
region from the nomograms using (0 O - 0j), (q 0 - q, ) and CFP. The mean surface sensible and 
latent heat fluxes may then be obtained by multiplying the mean sensible and latent heating by 
the mean wind speed in the boundary layer. 

The sensitivity of the solutions shown in the nomograms to the variations in the initial 
coastal soundings and large scale subsidence is also investigated. The general behavior of the 
mean sensible heating, the potential temperature and the height of the boundary layer as a func- 
tion of downwind distance within the cloud free region and their relations to several important 
parameters are analyzed in t*,ie Appendix. 

3. DESCRIPTIONS OF THE BOUNDARY LAYER MODEL 

The Stage and Businger (1981a, b) model was developed for the growth and evolution of 
the marine boundary layer during cold air outbreaks. This is a Lagrangian model. The equivalent 
potential temperature and the total water mixing ratio are assumed to be well mixed in the 
boundary layer and are predicted from conservation equations including fluxes from the sea sur- 
face and entrainment across the top of the boundary layer. The entrainment rate is determined 
from the turbulent kinetic energy budget. When clouds form, infrared cooling at cloud top and 
heating at cloud base are also included. The reality of the model was tested against data gathered 
during the International Field Year for the Great Lakes (IFYGL). Since our main interest in 
this paper is the air-mass modification in the cloud free region, only the significant physical pro- 
cesses in the same region are briefly discussed here. The reader is referred to Stage and Businger 
(1981a, b) for details of the model. 

The equations for the cloud free boundary layer are given by 

dfl/dt = (PPP |„ + W t A0)/Z B (I) 
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( 2 ) 


ilq/clt * (w' q' l 0 + W e Aq)/Z B 

dZ„/dt * w e + W, (3) 

where 0, q, and Z B are the potential temperature, specific humidity and height of the boundary 
layer, respectively, d( )/dt the rate of change following a column of air moving at the mean vel- 
ocity of the mixed layer, w 'O' | 0 and w'q' ! 0 the surface turbulent fluxes, W c the entrainment 
rate, Wjthe lifting rate of the inversion base due to large scale convergence, and A0 and Aq cor- 
respond to the jumps in 0 and q across the base of the inversion. Note that Wj « -DZ B , where 
D is the large scale divergence and is a prescribed parameter in the model. We shall see that W, 
generally appears to have little effect on our results. However, the entrainment rate, W e , can 
have considerable influence on the boundary layer through the introduction of warm, dry air 
from the inversion to the boundary layer. Therefore, a proper determination of W c is crucial in 
any mixed layer model. Here W e is determined from the turbulent kinetic energy budget; i.e,, 
the energy source for the entrainment, which works against buoyancy, is derived from the turbu- 
lent kinetic energy of the boundary layer, and this in turn is generated primarily by buoyancy 
production. In the cloud free region, the surface sensible heat flux plays a very important role 
in determining W e . 

The surface turbulent fluxes are parameterized according to the bulk transfer formula as 

w 7 ? 7 | 0 « C T U(0 o -0) (4) 

w'q' | Q = C q U (q 0 - q) (5) 

where 0 Q and q 0 are the sea surface temperature (SST) and the corresponding saturation mixing 
ratio, respectively, U the mean velocity of the mixed layer, and C T and C q the bulk transfer co- 
efficients for heat and moisture, respectively. C q is taken to be equal to C T . For the parameter- 
ization of the unstable surface layer, the transfer coefficients are independent of wind speed and 
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dependent only on the virtual potential temperature difference between the sea surface and the 
mixed layer. Tins is because the transfer coefficients depend only on the ratio of surface rough- 
ness to surfuce layer height, out both are proportional to wind stress. This has the helpful con- 
sequence that the downwind profiles of 0, q, and Z B arc independent of wind speed out to the 
edge of the cloud (Stage, 1979); the wind speed determines only the travel time for the required 
ftlr-mass transformations. Although some authors have shown transfer coefficients to vary with 
wind speed, the variation is only slight for unstable conditions such as those under consideration 
(Kondo, 1975). 


4. SENSIBLE AND LATENT HEATING NOMOGRAMS 

The mean sensible and latent heating between the shore and some downwind distance in 
the cloud free region are defined as 


x 


S - pep x" 1 

/ 

0 

X 

Cj (Aq — 0) dx 

(6) 

E = pL x" 1 

/ 

0 

c q (<1 0 “ 9) dx 

(7) 


where x is the distance from the shore along the surface wind direction, p the air density, c p the 
heat capacity of air at constant pressure, and L the latent heat of condensation. It can be seen 
from (6) and (7) that S and E are the mean heating of the air column by the warm sea per unit 
distance of travel as a result of the surface sensible and latent heat fluxes, respectively. It can 
also be seen from (4) - (7) that the mean surface sensible and latent heat fluxes are equal to SU 
and EU, respectively. Combining (1), (2), (6) and (7), S and E can be expressed as 

S = fiCf x-' JT (1 + A,,)- 1 Z„ d& (8) 

0 

£ » pi x-' / (1 - A.)’ 1 Z„ dq (9) 

o 


where 

A e = W C A d^O' | 0 


( 10 ) 
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(II) 


Aq * -W c Aq/w'q' | y 

Here the entrainment coefficients A q and A q are the ratios of the sensible and latent heat fluxes, 
respectively, at the top of the boundary layer to those at the surface. Without these cross- 
inversion fluxes, tlie heat und moisture budgets of the boundary layer would be in gross error. 

As noted earlier wind speed does not affect the CFP or the downwind profiles of 0, q, and Zjj 
in the cloud free region. It can be seen from (8) - (1 1) that S, E, A#, A q , and the nomograms 
of sensible and latent heating derived with the associated model are independent of wind speed. 

In order to construct the nomograms for S and E, several numerical experiments correspond- 
ing to the initial soundings shown in Table I and 0 O « 6, 8, ..... . 2b°C for nondivergent flows 

were performed using the Stage and Businger model. The initial soundings were mainly taken 
from Stage. The relative humidities (RH) shown correspond to those at the surface. If there is 
an SST gradient, S and E can be determined from the nomograms using the mean SST oi the 
region of interest for 0q. 

Fig. 2 shows the mean column sensible heating (S) between the shore and 100 km (within 
the cloud free region) versus the land air-sea surface temperature difference (0 Q - 0j) for non- 
divergent flow with the initial surface relative humidity of 20% and 0 j * 3°C. The initial poten- 
tial temperature lapse rate (fy) is 3.8°C km" 1 . It can be seen that S is nearly linearly propor- 
tional to (0 O -0,). This suggests a new kind of parameterization, i.c., the mean sensible heating 
is proportional to the horizontal difference rather than the vertical difference of temperature. 

The slight nonlinearity is mainly due to the weak dependence of C T on (0 O - 0f); the change in 
C T is about 2% per degree change in (0 O - 0,). This result, which was evidently overlooked by 
Stage (1979) and Stage and Businger (1981a, b) because they were not concerned with heat 
fluxes, has important ramifications in measuring and modeling coastal air-sea interactions since it 
depends on measurable parameters (0 O - 0j) rather than on vertical profiles which are often not 
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available. On the other hand, since this finding is itself a result of a model parameterization, it 
needs to be verified experimentally, 

Similarly, Fig, 3 shows the mean column latent heuting (E) between shore and 80 km (with- 
in the cloud free region) versus the land air-sea surface humidity difference (q 0 »• qj) for non- 
divergent Hows, The solid line corresponds to (0 q - 0j) * I l°C, and the dashed line to (0q -0 j) 

* 7 °C, both with r 0 * 3,8°C km" 1 , It can be seen that E is linearly proportional to (q 0 ~ q j ) 
and slightly dependent on (0 O - 0j) through the C q (* C-j) dependence, This result also suggests 
a new parameterization in that the mean latent heating is proportional to the horizontal differ- 
ence rather than the vertical difference of vapor mixing ratio. Similar comments apply here as 
were made in connection with the temperature parameterization. 

Following the ideas of Figs. 2 and 3 and the analysis shown in the Appendix, the nomograms 
for estimating the mean sensible (S) and latent heating (E) between shore and the downwind 
cloud free distance (x) under nondivergent conditions are given in Figs, 4 and 5 (solid curves). 

The solid curves of Fig. 4 show the relation of S to (0 O - 0j ) and x, while those of Fig. 5 show 
the relation of E to (q 0 ~ qj) and x with weak dependence on (0 O - 0j ), The application of 
Figs. 4 and S is straight-forward. To obtain the column sensible heating for the cloud free region 
one enters Fig. 4 with the cloud free path (x = CFP), and the difference between the mean sea 
surface temperature (0 O ) and the initial shore surface air temperature (0j). The SST is readily 
obtained from the 1 1 Mm channel of the AVHRR on TIROS-N or the same channel on a geo- 
synchronous satellites. In the usually dry cold air outbreaks, little correction is necessary for in- 
tervening water vapor. If there is a gradient of SST along the wind direction one may use the 
average SST between shore and the cloud edge since the heating is essentially linear with (0 O - 0j)* 
The value of 0 , may be obtained from conventional measurements at the shore or from IR mea- 
surements of the ground skin temperature empirically corrected to shelter level air temperature. 
Since the nomograms are produced with a specific value of r^(- 3.8°C km" 1 ), a correction is 
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also required Ibr other temperature lapse rates at the shore, which should be obtainable from 
conventional measurements. This correction Is discussed in Section 5. Note that Figf, 4 and 5 
can also be applied to any region with x < CFP provided that 0 O Is the mean SST for the region 
concerned. 

The use of Fig. 5 for latent heat is analogous to that of Fig. 4. Here qg is the saturation 
mixing ratio at the SST of 0 Q . Care must be exerted here in the event of a gradient in SST with- 
in x because of the non-linear relation between q 0 ?.nd 0 Qt The humidity of the initial shore 
surface air (qj) can be obtained either from direct measurements or from the daslie I curves of 
Fig. 4, which is to be described shortly. As noted earlier. E depends slightly on (Oq - 0j) 
through the C q dependence (Fig. 3). Therefore, a correction is required for E to account for 
the (0 Q - 0j ) dependence. This correction factor may be written as 

E * E n (1 + 0.02 60) (12) 

«» - -Vn < l3 > 

where the subscripts N and obs denote, respectively, the nomogram and observed values. The 
heat fluxes for the cloud free region may be obtained simply by multiplying the column heating 
by the average mixed layer wind speed in the cloud free region. The latter may be estimated by 
near surface measurements at the shore. 

The dashed curves of Figs. 4 and 5 correspond to the CFP for various relative humidities of 
the initial shore surface air temperature of 3°C. Although large scale subsidence and the initial 
lapse rates of temperature and humidity at the shore also affect the CFP, the dashed curves of 
Fig. 4 suggest that, to a first approximation, CFP and (0 o - 0j ) can be used to determine qj 
for a known 0 j. As mentioned above q, is required in using Fig. 5 to obtain the latent heating 
if one does not have a direct humidity measurement. Indeed, because of the notable errors in 
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radiosonde measurements of humidity, this approach may turn out to be more accurate once the 
theoretical results are validated experimentally, 

Note also that the cloud free path should be measured by a geosynchronous satellite either 
in the visible or !R because it varies diurnally with the land air temperature (see Section 6). Both 
Figs. 4 and 5 are valid only fo i x within the cloud free region. While Fig, 4 shows that sensible 
heating is not especially sensitive to x (the variation in S is about 10% per 100 km uncertainty 
in x for a constant (0 Q - 0 , ) curve), one needs the CFP to determine the average SST (0 q) in 
the cloud free region when there is an i?ST gradient. 

An examination of Figs, 4 and 5 shows the following key features within the cloud free 
region: 

(1) For constant x, S increases essentially linearly with increasing (0 q - 0j) and E with 

(% 

(2) For constant (0 O - 0j). S decreases by 30 to 35% for x ranging from 0 to 300 km. 

For constant (0 O - 0j) and (q 0 - q j ), E also decreases with increasing x, As the dry cold air 
travels over a warmer sea with a constant SST, the air becomes warmer and wetter and the local 
air*sea temperature and humidity differences are gradually reduced; hence S and E decrease 
gradually. 

(3) For the same initial temperature and humidity soundings, an increase in SST causes an 
increase in S and E but a decrease in CFP; i.e., fluxes increase and clouds form sooner. 

(4) For constant 0j and 0 Q , a decrease in the relative humidity of the initial land air causes 
an increase in CFP and E of the cloud free region but a decrease in S of this region. 

Although the nomograms arc generated with 0j = 3°C, its generality can be seen from the 
analysis shown in the Appendix and the (numerical solutions shown in Fig. 6, The solid curves of 
Fig. 6 show the downwind variation of the mean sensible heating, the local vertical air-sea tem- 
perature difference, and boundary layer height for the case of 20% RH with 0 j = 3°C, (0 Q - 0j) 
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« 1 1°(’ mu! V 0 * d,8°C km" 1 . For this ease the cloud free path is 280 km. The symbols along 
the curves show values of S, (0 O ** 0), and at the cloud euge for other initial temperature and 
humidity conditions of the land air but with the same values of (Oq - 0|) and F^, We note that 
the initial temperature soundings for the 80% and 60% Kil cases and the initial humidity sound- 
ings for the cases of 0 { * -10° and 0 j * -7°C are the same as those for the case of 20% HU and 
0 , * 3 °(\ 


It can be seen that the profiles of (Oq - 0), S and Z B in the cloud free region are virtually 
independent of the initial temperature and humidity soundings provided that (0 q ~ 0| ) and Vq 
are fixed. (The effects of Vq are discussed in Section 5). However, the cloud free path indicated 
by the position of the various symbols is highly sensitive to the initial soundings, i,e„ clouds form 
sooner at higher Rll, This confirms our previous findings, namely Unit the mean sensible heating 
for the cloud free region is uniquely determined by (Oq - 0,), Vq, and OFF. We shall see that 
this is also consistent with the analysis in the Appendix. The fact that the boundary layer growth 
in the cloud free region is almost totally controlled by the surface sensible beat flux, and thus 
also by (Oq - 0 | ) and Vq , is also significant since the boundary layer height is a parameter which 
should also be measurable from space. Indeed, other things being equal, the growth rate of the 
boundary layer and sensible heating may be used ns proxies for one another. Fq» (AS 5) in the 
Appendix shows fl at the sensible heating is directly proportional to the average value of d(7. B -)/ 
dx in the cloud free region. We shall discuss this further below. 


Our analyses also indicated that the mean latent heating for the cloud five region may be 
determined nearly uniquely from a knowledge of (q 0 -q, ), (0q - 0 |), Vq and CFP. For a fixed 
(0 q - 0,), (q 0 - q, ) and Vq, the mean latent heating in the cloud five region was found to be 
slightly dependent on the initial humidity lapse rate, 1'. However, the effect is very small, as 
can he seen from Table 2, which shows the percentage errors in S and E at 100 km from the 
coast due to errors in various factors. The percentage errors in S and E due to an error of 1°C 


in (0 O - 0| ) and of I g kg" 1 in (q 0 - q, ) arc smaller for larger (0 O - 0j ) and (q 0 - q j ), respective- 
ly. They amount to about 1 1% in S for (0 Q - 0 j ) = 9 - II °C and 15% in E for (q 0 - qj ) s 6 - 
8 g kg”* . A 10 km uncertainty in CFP causes about a 1% and a 0.5% error in S and E, respec- 
tively. A I0" 5 s" 1 uncertainty in the divergence causes about a 0.2% and a 0.03% error in S and 
E, respectively. This suggests that the nomograms (Figs. 4 and 5) generated from the nondiver- 
gent flow may generally be used to estimate S and E for divergent flows. 

5. EFFECTS OF INITIAL POTENTIAL TEMPERATURE LAPSE RATE 

For a constant 0 j , an initial Ifuid air column associated with a large is warmer than that 
with a small ty, After modification from below by warmer water the former tends to maintain 
a warmer temperature in the boundary layer and thus produces a smaller driving force for the 
sensible heat flux. In spite of the smaller sensible heating with larger T# the boundary layer al- 
ways remains warmer because of its shallower depth as compared to that with small Vq , There- 
fore, for a fixed (0 O - 0j ), S decreases with increasing Similarly, for a fixed (0 Q - 0| ) and 
(q 0 — q | )» E also decreases with increasing 

For r 0 greater than 3.8°C km -1 , the sensible heating nomogram (Fig, 4) overestimates the 
sensible heating and thus requires a negative correction factor, and conversely. According to our 
analysis, the Z B -normalized correction factor for the mean sensible heating may be approximated 
as (see Appendix) 

(AS/S)/Z B - -0.25 Id + A e )/(l + 2A 0 )j [(r 0 - 3.8)/(0 Q - 0j>] (14) 

Here Z B is the boundary layer height at the downwind edge of the region in which S is computed, 
and Aq the entrainment coefficient of potential temperature defined in (10) as the ratio of sensi- 
ble heat flux at the top of the boundary layer to that at the surface. The value of Aq computed 
from Stage and Businger’s model is generally about 0.3 in the cloud free region for (0 Q -0j) > 
7°C. Note that the entrainment coefficient of virtual potential temperature is equal to 0.2 in the 

13 


O 


same model. These two entrainment coefficients are in reasonable agreement with those obtained 
in other studies (e.g., Tennekes, 1973; Willis and Dcardorff, 1974; Randall, 1980), Figs, 7 and 8 
show the ^-normalized correction factors for the mean sensible heating, (AS/S)/Z B , and the 
mean latent heating, (AE/E)/Z B versus F 0 for various (0 Q - 0j) as derived from the results of 
numerical experiments, It can be seen that Fig. 7 is consistent with (14) and that the Z B -normal- 
ized errors increase with increasing deviation of from 3.8°C km" 1 , and decrease with increas- 
ing (0 O - 0j ). The correction factors (AS/S) and (AE/E) for the cloud free region may be ob- 
tained by multiplying the ordinates by Z B (km) at the cloud edge. The boundary layer height 
may be measured by a variety of ground based remote sensors such as lid ar or sensitive micro- 
wave radars, to the future, it should be measurable by a spaeeborue Ildar scattering off the sea 
salt aerosols at the top of the convective boundary layer, 

According to our analysis, the boundary layer slope within the cloud free region under non- 
divergent conditions may be written as (see Appendix) 

dZr,/dx * E(I + 2A 0 )(0 0 - 0, )/( V Q Z B ) - d + A#)] C x (15) 

For large (Oq ~ 0j ) and small F^ and Z B , (15) may be approximated as 

dZ B 2 /dx ~ 2(1 + A^)C T (0 q - 0, ) IV 0 (16) 

Eq. (16) suggests that Vq may be estimated from a knowledge of d(Z B 2 )/dx and (0 Q - 0 , ), Fig, 

9 shows the slope of Z B 2 versus for various (0 q -*0|) under noudivergent conditions as de- 
rived from the results of the numerical experiments. Here d(Z B 2 )/dx is determined from the 
first 50 km within the cloud free region (or from the cloud free region if the CFP is smaller than 
50 km). It can be seen that d(Z B 2 )/dx increases with increasing (0 Q - 0j ) and decreases with 
increasing Vq , which is consistent with (16). For nondivergent flow, the potential temperature 
lapse rate for the unmodified air (i,e„ above the convective boundary layer) at time t, 7 ^(t), is 
the same as that of the initial sounding, P#, For divergent How, 7^(0 is gradually increased by 


1 


the divergence and may be given as (Carson, 1973) 

70 <0 * r 0 exMDt) (17) 

where D is the divergence. Therefore, divergence tends to decrease Z B by increasing the poten- 
tial temperature lapse rate. However, the effects of divergence on d(Z B 2 )/dx are expected to be 
very small for the first 50 km or the first hour of travel. 

6. OBSERVATIONS 

Fig. 10 shows a T1ROS-N AVHRR 1R image for a cold air outbreak bn 17 February 1979 
at 1941 GMT. The cloud streets are aligned essentially along the wind direction, from the north 
and northwest, as seen from the winds at the three buoys as sown in Fig, 1 1. Initial cloud for- 
mation typically occurs about 100 km off the nearest upwind shore so that the boundary of the 
clouds is virtually an image of the eoastline, South of Long Island, New York, the image shows 
surface water temperatures increasing from about -1,5°C at the shore to 2°C about 50 km south. 
Approximately 500 km south of Long Island (36°N latitude), the clouds become solid and are 
extremely uniform in the east-west direction over a distance of about 1 100 km. Note that the 
coldest cloud top temperatures of -21°C (dark green) occur at latitude 38°N where the clouds 
remain in the form of streets, Although it is almost certain that the cloud tops continue to grow 
southward of the coldest point, the top temperatures grow progressively warmer. This is due to 
the growth of the clouds into the inversion layer so that higher clouds are associated with warmer 
temperatures. 

Fig. 1 1 shows the surface weather map at 1500 GMT (1000 EST), about 5 hours earlier 
than the time of TIROS-N AVHRR picture (Fig. 10). A high pressure system centered in the 
Great Lakes area dominates the entire east coast of North America. Cloud street orientation is 
generally observed to be along the mean wind in the planetary boundary layer and about 1 5° to 
the left of the surface gecstrophic wind (Lemonc, 1973). This is also consistent with the case of 
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17 February 1979, Fig. 12 shows the weekly mean SST for the period of 14-21 February 1979. 
The temperatures are seen to be slightly warmer than the AVHRR 1R brightness temperatures 
for the cloud free coastal region, Note that the air flows across the SST isotherms for the region 
from Maine to Long Island, where the cloud streets arc best organized. Further south, the air 
flow is nearly parallel to the SST isotherms. Also note the moderate snowfall (Fig. 1 1) reported 
at Cape Hatteras, the snow quite obviously having originated in the low level clouds formed over 
the relatively warm sea. 

In order to examine the present method the case of 17 February 1979 was fully analyzed 
and the sensible heating derived from the nomogram (Fig. 4) was then compared to that com- 
puted from the bulk transfer formula, A buoy located at about 75 km downwind from the Long 
Island shore (40.1°N, 73°W) and several coastal stations provided hourlj air (T a ) and sea surface 
temperatures (T s ) and wind speed (U) for the computation of sensible heat flux (SH) with the 
bulk transport formula, which is 

SH = pc p C H (T s - T U )U (18) 

Here Cj^ is the bulk transfer coefficient for the sensible heat flux, and is taken to be 1.5 x 10“ 3 
(Rondo, 1975) for the diurnal ranges of temperature and wind in this case. Note that this value 
of C H is consistent with that of C T generated in Stage and Businger’s model, which is about 1.5 
'v 1,6 x 1CT 3 for this kind of situation ((0 O ~9j) ^ 15 - 19°C). The sensible heat fluxes at the 
shore and at the buoy were computed from (18). The mean sensible heating was then computed 
by dividing the mean sensible heat flux by the mean wind speed between the shore and the buoy. 

The cloud free path, averaged over about 120 km wide coastal region south of Long Island, 
was measured manually from the hourly GOES-2 IR cloud pictures. Because the mean sensible 
heating for the cloud free region discussed in Section 4 is the average along the trajectory in a 
Langrangian sense, and the manually measured CFP is fairly subjective, a three-point running mean 


was applied to all the data (CFP, U, T s and T a ). In order to estimate the mean sensible heating 
for the cloud free region from the nomogram, we need both (0 O -0j) and the CFP, Since it 
takes time for the initial land air to form the first cloud downstream, we used values of 0 1 ob- 
served at the shore at times earlier than the cloud observations differing by the travel time from 
shore to cloud edge corresponding to the mean wind speed between the shore and the buoy. The 
mean SST for the cloud free region, 0 O , was determined from the CFP and the constat SST field, 
which was estimated from the SST at the coastal station and the buoy. The coastal SST field 
was found to have a gradient of about 0.07°C knf 1 and was rather steady. 

Fig. 13 shows the diurnal variations of cloud free path (circles) and the associated mean sen- 
sible heating for the 17 February 1979 cold air outbreak south of the Long Island, New York 
coast. The solid curve is the mean sensible heating between the shore and the buoy (75 km) 
computed from the bulk transfer formula, while the squares represent the mean sensible heating 
between the shore and the cloud edge (60 to 120 km) estimated from the sensible heating nomo- 
gram (i.e., Fig. 4). Note that the errors in the nomogram estimates due to divergence and the 
initial potential temperature lapse rate (fy) are probably fairly small because of the small diver- 
gence (1.5 x 10" s sec -1 ) and large values of (0 q - 0 ( ) (about 15 - 19°C for most of the data) 
and the small deviation of Vq from 3.8° km" 1 (r# is about 1.3 - l .8°C km" 1 ), 

We see that the results obtained from the present method are generally in very good agree- 
ment with those from the bulk transfer formula. Since there is no diurnal variation in the coast- 
al SST field as evidenced both at the shore station (Fire Island) and the buoy, the variation in 
the mean sensible heating between the shore and the buoy is due mainly to the variation in the 
initial shore surface air temperature (0|). The initial shore surface air temperature varies from a 
minimum of -18°C around 0700 EST to a maximum of -I0°C near 1500 EST. Note that the 
maximum of S (between shore and the buoy) generally corresponds to the coldest land air tem- 
perature and vice versa. On the other hand, there is a slight diurnal variation in the mean SST 


17 


over the cloud free region (0 O varies front 1,2 to 3.8°Q because of the gradient of SST (about 
0.07°C km" 1 ) and the variation in CFP, This is due only to the fact that the varying cloud edge 
is the outer limit over which 0 Q is averaged. In short, the variation in S in the cloud free region 
is due to the variation in tiie initial shore surface air temperature, 0 ( , and the change in the ef- 
fective SST, 0 O . For a CFP smaller than 75 km, S of the cloud free region (squares) is generally 
smaller than that of the region between the shore and the buoy (solid curve), because the effec- 
tive 0 O is smaller in the former than in the latter. For CFP larger than 75 km, the reverse is 
true. This accounts for most of the difference between the present and the conventional methods 
of estimating heat flux. 

The daily averages of S and SH computed from the bulk transfer formula are about 25 J in" 3 
and 249 W nf 2 , respectively, with the observed daily mean wind speed of aoout 10 ms” 1 . These 
values agree almost identically with the average values obtained from the nomogram. Although 
Fig, 13 indicates that S is overestimated after 1200 GMT for the reason noted above, these were 
balanced by the underestimates prior to that time, thus giving excellent agreement for the daily 
averages. 

The present method was also applied to the case of the cold air outbreak analyzed by Henry 
and Thompson (1976). The values of 0 O , 0, , q 0 , q, , CFP, and U were chosen to be 20°C, -4°C, 
14.5 g kg -1 , 1 gkg' 1 , 100 km and 10 ms" 1 , respectively. The corresponding sensible and latent 
heat fluxes estimated from the nomograms were 400 and 600 W m" 2 , respectively, for the cloud 
free region with a travel time of about 3 hours. These values are much smaller than those ob- 
tained by Henry and Thompson. The sensible and latent heat fluxes are about 27% and 72%, 
respectively, of their corresponding values. Since the temperature and humidity soundings near 
the cloud edge used by Henry and Thompson were estimated, not observed, it is likely that the 
estimated soundings were too warm and too wet, thus producing excessive apparent sensible and 
latent heat fluxes. Also, the heating and drying effects of entrainment through the top of the 


boundary layer, which are important in the heat and moisture budgets, are not included in the 
Henry and Thompson budget calculation. This is surely another m^jor cause of the discrepancy. 

1 , CONCLUSIONS AND DISCUSSION 

In the present paper we have shown how the sensible and latent heating of the atmosphere 
due to surface fluxes is related to the downwind distance offshore at which clouds first form 
during cold air outbreaks over the warmer coastal waters. This distance, the so-called cloud free 
path or CFP, and the mean sea surface temperature, both of which arc measurable from space, 
along with the temperature and moisture content of the initial shore surface air, both measurable 
by conventional means, are all that are required to utilize the nomograms (Figs. 4 and S) devel- 
oped here. Other things being equal, the greater the land air-sea surface temperature contrast, 
the greater the heat and moisture fluxes and the closer to shore the clouds will form. 

Tlic results are not sensitive to flow divergence or to initial lapse rate of humidity at the 
shore, but are dependent upon the initial lapse rate of potential temperature, fy, If this is not 
known from conventional soundings at the shore, it may be deduced from the downwind slope 
of the square of the boundary layer height (Z B 2 ) which is proportional to the land air-sea sur- 
face temperature difference and inversely proportional to Fg. The mean slope of (Z B 2 ) is also 
directly proportional to sensible heating so that it may be used as a proxy for the latter. The 
slope of the boundary layer height may be measured by a variety of ground based remote sen- 
sors such as lidar, acoustic radar, or sensitive microwave radars, In the future, we expect that it 
will also be measurable by spaceborne lidar scattering off the sea salt aerosols at the top of the 
convective boundary layer. Another important future step is to extend the present approach fur- 
ther from the coast beyond the cloud edge because it is clear that significant heating continues 
to occur after clouds have formed (Henry and Thompson, 1976). 

Within the cloud free region, the profiles of potential temperature, water vapor mixing ratio 
and the height of the boundary layer along the wind direction were found to be independent of 


wind speed. This also makes the downwind profiles of the mean sensible and latent heating in- 
dependent of wind speed, The generality of these results needs to be investigated further through 
the comparison with other models (e.g„ Brost, 1976, Randall, 1980, Moeng and Arakawa, 1980) 
and empirical studies, The various models may have different values of the entrainment coeffic- 
ients, Theoretical analyses (see Appendix) suggest that the growth and heating of the boundary 
layer within the cloud free region are not sensitive to the variation in the entrainment coefficient 
(A q) so long as it is within a reasonable range, If the parameterization of the turbulent transfer 
coefficients through the air-sea interface is similar to that in Stage and Busiuger’s model, it is 
expected that the model dependency of the nomograms will be very small. Nevertheless, we rec- 
ommend validation experiments and fine tuning of our results before they are applied more 
broadly, 

The method was tested in two cases, The first one showed good agreement in sensible heat- 
ing as compared with that computed from buoy measurements and the bulk transfer formula. 
However, a major discrepancy was found in our comparison to the case reported by Henry and 
Thompson (1976) due to factors discussed earlier. 

There arc some interesting by-products of this work which are worth noting. First of all, 
the model results indicate that the sensible heat flux is almost linear with the land air-sea surface 
temperature difference, and that the latent heat flux is linear with the land air-sea surface vapor 
mixing ratio difference, the latter with a small correction for the temperature difference. This 
suggests the possibility of a new kind of flux parameterization for the convectively unstable case 
which depends upon the horizontal differences across the coast instead of upon the vertical dif- 
ferences as in the classical bulk transfer equations. 

Another interesting aspect of the present work is that it appears that the method can be 
used in reverse to check the validity of some of the basic assumptions and features of boundary 
layer models. For example, it would be useful to determine whether or not the CFP and the 
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other parameters ( 0 , q, Z B , S and E) on which it depends are indeed independent of the wind 
speed as the Stage and Businger (1981a, b) model suggests. In addition, observations such as 
these should permit an indirect determination of the overall validity of models such as that used 
here, While it may not be possible to separate the effects of various processes or parameters in 
the boundary layer model, the position of the first cloud formation should provide fairly sensi- 
tive tests of its overall performance. 

Finally, it is evident that the physical processes and results discussed here are related to the 
lake-effect snowstorms observed on the leeward side of lakes during cold outbreaks (Lavoie, 1972; 
Passarelli and Braham, 1981). It seems likely that our method should be adaptable to the now- 
casting of such snowfalls since the snowfall rate is related to the net moisture flux, or the differ- 
ence between that coming from the lake (ocean) surface and that transported across the inversion. 
It is of interest that both these important moisture supply and loss terms are automatically in- 
corporated in determining the cloud free path, 
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APPENDIX 

Analytical Formulation for the Cloud Free Region 

The heat budget equation lor the cloud tree boundary layer may be written as 

d0/dt * | 0 4- W e A0)/Z B (Ai) 

where 0 and Z B are the potential temperature and height of the boundary layer, respectively, t 
time, w'0' | 0 the surface turbulent flux of potential temperature, W e the entrainment rate, and 
A 0 the jump of 0 across the inversion base, For nondivergent conditions (dZ B /dt * W„), (Al) 
may be rewritten as 

d0/dZ B * f (1 + A 0 )/A 0 J (A0/Z b ) (A2) 

where the entrainment coefficient is given by 

A 0 = W c A Q/W' | 0 (A3) 

The jump in potential temperature at the inversion base may be given as 

A0 “ 0j + Tfl Z B - 0 (A4) 

where 0 j is the initial surface potential temperature at the shoie and T# the initial potential 
temperature lapse rate. Using (A4), (A2) may be rewritten as 

d0/dZ B = [ (1 + A e )/A 0 ) (F 0 - (0 - 0, )/Z B J (A5) 

which may be integrated to give 

0 - 0(O)M + 0, (I - M) + [ (1 + A 0 )l(\ + 2A 0 ) ] V d IZ B - Z B (0) M ] (A6) 

M = [Z B /Z B (0)1 - [ < i+ V /a 0 i 
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where 0(0) and Z B (0) are the initial potential temperature and height of the mixed layer, respec- 
tively, Note that the initial temperature sounding is generally stable and Z B (0) is very small, For 
the convectively unstable boundary layer, Aq is generally about 0.3 in the cloud free region and 
Z B is generally much larger than Z B (0). Therefore (A6) becomes 

0 * 0, + 1(1 + A 0)1(1 + 2A e )\ r 0 Z B (A7) 

Eq, (A7) is a linear relationship between the potential temperature and height of the convective 
boundary layer in the cloud free region, 

Eq, (Al) may also be written as 

d0/dt * (1 +A0)W 7 0 7 | O /Z b (A8) 

where 

wV to * C T U(0 O - 0) (A9) 

C T is the bulk transfer coefficient, U is the mean wind speed in the boundary layer, and 0 O is 
the sea surface temperature. 

Combining (A8) and (A9), results in 

d0/dx « C T (1 + Aq) (0 o - 0)/Z B (A10) 

Where x is the distance from the shore along the wind direction. 

Combining (A7) and (A10), results in 

dZ B /dx = 1(1 +2A 0 )(0 O -0j)/(r^Z B ) - (I + A^)l C T (All) 

For large (0 Q - 0j ) or small and Z B , (All) may be approximated as 

dZ B /dx = C x (l + 2A 0 ) (0 O - 0i)/(r 0 Z B ) <A12) 
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which may be integrated to give 


Z B : " z„ 2 (0) + 2 C T (I + 2A 0 ) (9 0 - 0, ) x/l'o (A13) 

where 

x 

Op * x ^ / Cj dx 

0 

The mean sensible heating between the shore and some downwind distance in the cloud free re- 
gion is defined as 

x 

S * pep x _1 / C T (0 o ~ 0)dx (A 14) 

0 

Using (A7) and (A 10), (A 1 4) may be expressed as 

. x 

S = pep x” / [0.5r 0 /(l + 2A 0 )J fdz B 2 /dx] dx (A15) 

0 

Eq. (A15) indicates that the sensible heating in the cloud free region is directly proportional to 
the average value of d(Z B 2 )/dx in the same region, This supports the earlier statement that the 
sensible heating and the growth rate of the boundary layer arc proxies for one another, Using 
(All), (A15) may also be expressed as 

S = p Cp C T (0 o -O I )-O.5pc p C T r 0 Z B (l + A 0 )/ (1 4- 2A 0 ) (A16) 

From (All) and (A16) we see that S can be uniquely determined by (0 Q - 0 j), V 0 and x, Note 
also that the second term of the R.H.S. of (A16) is generally much smaller than the first term, 
Using (A 13) and (A16). the fractional error in S due to the change in r 0 rnay be approximated as 

(dS/dr 0 )/S ± "0.25 1(1 + A 0 )/(l + 2A 0 )J IZ B /(0 O - 0,)) (A17) 

As mentioned earlier, A# computed from Stage and Businger’s model is generally about 0.3 in 
the cloud free region for (0 O - 0j ) > 7°C, It can be seen from (A7), (A13) and (A16) that an 
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increase of A q from 0,2 to 0,4 causes a 13% increase in Z B , a 3% increase in (0 - 0 j ) and a 3% 
decrease in [S - pc P C T (0 O - 0, ». This suggests that our results are insensitive to A# and that 
the model dependency of the nomograms is almost completely controlled by the parameteriza- 
tions of the turbulent transfer through the air-sea interface, particularly upon the transfer coeffic- 
ients, 
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TABLE 1 . Initial soundings 
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TABLE 2. Errors in sensible and latent heatings at 100 km from the shore. 
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FIGURE CAPTIONS 


Figure I, Cold air outbreak on 3 January 1979 at. seen by GOES-2 visible channel. Note that 
shoreward edge of clouds is nearly an image of the coastline. 

Mean column sensible heating (S) between shore and 100 km versus land air-sea sur- 
face temperature difference ( 0 O - 0 j ) for an initial surface RH of 20% with 0 j * 

3°C and initial potential temperature lapse rate (fy) of 3.8°C km" 1 . 

Mean column latent heating (E) between shore and 80 km versus land air-sea surface 
humidity difference (q Q - qj ). Solid line corresponds to (0 Q - 0 j ) = 1 1°C, and 
dashed line to (0 O - 0 j ) = 7°C. All for the initial potential temperature lapse rate 
(r 0 ) of 3.8°C km” 1 . Symbols in legend correspond to various initial soundings. 

Mean column sensible heating (S) between shore and downwind cloud free distance 
(x) versus x for various land air-sea surface temperature differences (0 q - 0j) under 
nondivergent conditions and initial potential temperature lapse rate of 3.8°C km" 1 
(solid lines). Dashed lines correspond to the cloud free paths for various surface 
relative humidities (RH) and initial shore surface air temperature (0 j ) of 3°C. Each 
cross indicates a numerical experiment. 

Figure 5. Sames as Fig. 4, except for the dependence of mean column latent heating (E) on 
land air-sea surface humidity difference (q Q - qj ) and land air-sea surface tempera- 
ture difference (0 Q - 0 j ). Equation at bottom shows corrections in E due to values 
of (0 O - 0j ) not plotted on nomogram (see text). 

Figure 6. Mean sensible heating (S) between shore and downwind cloud free distance (x), 

local vertical air-sea temperature difference (0 Q - 0), and boundary layer height (Z B ) 
versus x for 20% relative humidity (RH) with 0j = 3°C (solid curves as marked). 


Figure 2. 


Figure 3. 


Figure 4. 
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Figure 7. 

Figure 8. 
Figure 9. 
Figure 10. 


Figure 1 1 , 

Figure 12. 
Figure 13. 


Symbols correspond to other initial shore conditions as shown in legend. All for 
(0 O - 0, ) « 1 1°C, and = 3.8°C km‘ l . 

Zp-normalized correction factor for mean sensible heating (AS/S)/Z B versus the in- 
itial potential temperature lapse rate (fy) for various land air-sea surface temperature 
differences (0 O - 0j). Z B (km) is the boundary layer height at the downwind edge 
of the region in which S is computed. 

Sames as Fig. 7, except for Z B -normalized correction factor for mean latent heating, 
(AE/E)/Z b , 

Slope of Z B 2 versus the initial potential temperature lapse rate (r^) for various 
(0 o - 0j) under nondivergent conditions. 

Cold air outbreak on 17 February 1979 as seen by Tiros-N AVHRR IR channel at 
1941 GMT. Brightness temperatures are color coded according to color bar at bot- 
tom in steps of 1.8°C. Abscissa and ordinate are marked in longitude and latitude. 
Coastal waters in pink, red, and purple from Maine to Chesapeake Bay are in range 
of -1.5°C to 5.5°C, Northermost cloud edges appear stippled because clouds are 
broken. Cloud streets appear clearly in green south of Long Island and Connecticut, 
and in blue or green southeast of coast of Maine. Coldest cloud temperatures are 
green (-21°C) along 38° latitude. 

Surface weather map of 17 February 1979 at 1500 GMT. Note snow at Cape Hat- 
teras and winds at three buoys. 

Weekly mean sea surface temperature for the period 14-21 February 1979. 

Diurnal variations of cloud free path (circles) and the mean sensible heating between 
shore and cloud edge (solid curve and squares) for 17 February 1979 cold air out- 
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break off coastal area of Long Island, New York. Solid line is computed from the 
bulk transfer formula and the squares are estimated from the nomogram, both using 
data at shore and buoy. 




LAND-SEA SURFACE TEMPERATURE DIFFERENCE <0 O -0i) (°C) 

Figure 2, Mean column sensible heating (S) between shore and 100 km versus land 
air-sea surface temperature difference ( 0 Q - 0 j) for an initial surface RH 
of 20% with 6 j = 3°C and initial potential temperature lapse rate (Va) 
of 3.8°C km" 1 . 



COLUMN LATENT HEATING, E (J m -3 ) 



LAND-SEA SURFACE HUMIDITY DIFFERENCE 

(q 0 -q,l(g/kg) 


Figure 3. Mean column latent heating (E) between shore and 80 km versus land air- 
sea surface humidity difference (q 0 - qj ). Solid line corresponds to (0 O - 
0j) * 1 1°C, and dashed line to (0 O - 0j) * 7°C. All for the initial poten- 
tial temperature lapse rate (fy) of 3.8°C km"* . Symbols in legend corres- 
pond to various initial soundings. 
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COLUMN SENSIBLE HEATING, S(J m~ 3 ) 



DOWNWIND CLOUD FREE DISTANCE, X (km) 


Figure 4. Mean column sensible heating (S) between shore and downwind cloud 
free distance (x) versus x for various land air-sea surface temperature 
differences (0 Q - 0 } ) under nondivergent conditions and initial poten- 
tial temperature lapse rate of 3.8°C km -1 (solid lines). Dashed lines 
correspond to the cloud free paths for various surface relative humid- 
ities (RH) and initial shore surface air temperature (6 j ) of 3°C. Each 
cross indicates a numerical experiment. 
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:NT HEATING, E n (J m-3) 



DOWNWIND CLOUD FREE DISTANCE, X (km) 

Figure 5. Same as Fig. 4, except for the dependence of mean column latent heating (E) 
on land air-sea surface humidity difference (qQ - qj) and land air-sea surface 
temperature difference (0 Q - 0|), Equation at bottom shows corrections in E 
due to values of (0 Q - 0j ) not plotted on nomogram (see text). 
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(0„-8)(°C) OR S (J m' 3 ) 



DOWNWIND DISTANCE, X (km) 


Figure 6, Mean sensible heating (S) between shore and downwind cloud free distance 
(x), local vertical air-sea temperature difference (0 Q - 0), and boundary 
layer height (Z B ) versus x for 20% relative humidity (RH) with 0 j « 3°C 
(solid curves as marked). Symbols correspond to other initial shore condi- 
tions as shown in legend. All for (0 O - 0j) = 1 1°C, and Tq ~ 3.8°C km” 1 . 
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BOUNDARY LAYER HEIGHT, Z B (km) 





INITIAL POTENTIAL TEMPERATURE LAPSE RATE, la CC/kml 


Fieure 7 Zn-normated correction factor for mean sensible heating (AS/S)/Z B versus 
the initial potential temperature lapse rate (ty) for various land air-sea surface temper- 
ature differences (»„ - 9,). Z B (km) is the boundary layer heigllt at the downwind 
edge of the region in which S is computed. 





Figure 8. Same as Fig. 7, except for Zg -normalized correction factor for mean latent 
heating, <AE/E)/Zg. 
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INITIAL POTENTIAL TEMPERATURE LAPSE RATE, f e (°C/km ) 

Figure 9. Slope of Z B 2 versus the initial potential temperature lapse rate (r#) for 
various ( 0 Q -8 j) under nondivergent conditions. 
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Figure 10. Cold air outbreak on 17 February 1979 as seen by Tiros-N AVHRR IR 
channel at 1941 GMT. Brightness temperatures are color coded accord- 
ing to color ' '*■ at bottom in steps of 1 ,8°C. Abscissa and ordinate are 
marked in longitude and latitude. Coastal waters in pink, red, and pur- 
ple from Maine to Chesapeake Bay are in range of -1.5°C to 5.5°C. 
Northermost cloud edges appear stippled because clouds are broken. 
Cloud streets appear clearly in green south of Long Island and Connect- 
icut, and in blue or green southeast of coast of Maine. Coldest cloud 
temperatures are green (-21°C) along 38° latitude. (\ 


44 





1O0 200 


KILOMETERS 


V 



U.S. EAST COAST 


Figure 1 1. Surface weather map of 17 February 1979 at 1500 GMT. Note snow 
at Cape Hatteras and winds at three buoys. 
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Figure 12. Weekly mean sea surface temperature for the period 14-21 February 

1979. 
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Figure 13. Diurnal variations of cloud free path (circles) and the mean sensible heating 
between shore and cloud edge (solid curve and squares) for 17 February 
1979 cold air outbreak off coastal area of Long Island, New York. Solid 
line is computed from the bulk transfer formula and the squares are esti- 
mated from the nomogram, both using data at shore and buoy. 
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